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(57) Abstract 

A method for the active control of sound radiated in a duct (I) comprising sensing sound from an array (4) of sensors located on the 
inside surface of a flow duct and controlling an array (3) of secondary sources (loudspeakers) so as to minimize sound radiated in the far 
field. This is done by controlling the sound sources so as to minimize a cost function. Depending on whether total sound power or sound 
power towards the sideline is to be reduced, different cost functions are used by the method. 
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I 

METHOD AND APPARATUS FOR THE ACTIVE CONTROL OF SOUND RADIATED FROM 

FLOW DUCTS 

The invention relates to the active control to limit harmonic sound radiated towards the sidelines from 
ducts containing a subsonic, uniform flow. It is in particular applicable to limitation of noise radiated 
from circular ducts such a gas turbine intake. 

There is growing interest in applying active noise control to reduce the fan tones radiated from aircraft 
turbo fan intakes, prompted by increasingly strict legislation regulating noise levels in populated areas. 
The traditional approach of sound proofing is to line the air inlet with sound absorbing material. 
However with the increasing trend towards shorter inlet required for the higher efficiency of high bypass 
ration engines, there is less space for the sound proofing material to be located. 

An alternative approach in recent years has been to investigate active noise control. Patent applications 
WO 95/19075 and WO 94/08540 both describe active controllers for flow ducts having internally located 
sensors and sound sources. No details however are given of the control algorithm. Patent US 5355417 
discloses a configuration for the active control of aircraft engine inlet noise by including an array of 
circumferentially arranged sound sources mounted inside an inlet duct as well as an array of sensors 
arranged in a ring. Again the active control algorithm is not disclosed in any great detail. Moreover the 
results of this system show an increase in sound propagated towards the sidelines, which are the 
important regions for sound reduction. 

Using external sensors, the inventors have determined a clear relationship between reductions in the 
radiated far field acoustic pressure in the radiated far field (in the region well away from the duct) and 
corresponding reductions in the internal transmitted field inside the duct. This provided relationships 
that allowed the use of internal sensors and internal sources to reduce noise levels in the far field towards 
the sidelines in a controlled and determined fashion. In other words an obtainable quantiry at the duct 
wall has been determined which is a function of acoustic pressure, which has a robust and stable 
relationship to the far field acoustic pressure for the important range of radiation angles that contribute 
most to the annoyance of those living beneath the flight path. 

It is an objective of the invention to provide both an active noise control arrangement and a method to be 
implemented which reduces the fan tones radiated from the intake of turbofan engines in the sidelines. 

ArmrHino io the invention is provided a duct for fluid flow having means for active control of sound 
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radiated therefrom, said duct comprising sound sensors located on the inner surface of said duct and 
grouped together in one or more planes transverse with respect to the duct axis, and at least one 
secondary source whose operation is a function of sound received at said sound sensors characterised in 
that the axial spacing of said transverse planes is not more than 0.5>. min (l+M 2maN ) where /^ tn is the 
wavelength corresponding to the radiated tone frequency of interest and M zmax is the maximum Mach 
number of the free stream flow in the duct. 

By using sensors external to the engines to observe directly the far field radiated sound the inventors 
have determined a method to controls loudspeakers, or so called secondary sources, so as to minimise 
engine noise in the far field. 

The invention also provides a method for the active control of sound radiated from a fluid flow duct 
comprising: 

a) sensing sound from an array of sensors located on the inside surface of the duct and 

b) controlling an array of secondary sources located on the inside surface of the duct, so as to minimise 
sound radiated in the far field in a pre-set band of angles to the duct axis, characterised in that only sound 
propagated within a pre-set angular interval to the duct axis is effectively sensed and used to control the 
operation of said secondary sources. 

Further the inventors have determined a relationship between the internal and external fields which has 
been incorporated into a cost function which, upon minimisation, has the desired effect of producing 
sound pressure level reduction in chosen radiation angles. Manipulation of the in-duct acoustic pressure 
as a result of the observation can be used to minimise transmitted sound and therefore far field sound. 

More specifically it relates to a line array of in-duct wall mounted discrete sensor elements whose 
pressure signals can be processed to provide an estimate of transmitted sound field with propagation 
angle. This measurement has been shown to be closely related to the variation in the radiated field 
versus polar angle. The internal sound field at the sensor elements is used for controlling the radiated 
field in the important range of radiation angles towards the sidelines. The important aspect of the induct 
error sensing methods proposed here is the simplicity, of the algorithm used. 

Advantages of the method are the simplicity and generality to all circular flow ducts since the proposed 
algorithms only require the fan tone frequency of interest and the speed of the free stream flow as mpui 
variables. 
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The invention will be described with reference to the following figures of which: 

Figure I shows a schematic diagram of a circular duct containing uniform axial flow, comprising an 
array of sensors and secondary sources. 

Figure lb shows the relationship between resultant wavenumber km. normal to the localwavefront. the 
propagation angle 8 mn and axial wavenumber k^ . 

Figure 2 shows the a typical variation of axial wavenumber with propagation angle with the ideal 
receiver response. 

Figure 3 shows the directivity function of a ten element liiie array at the design frequency steered at 45° 
for zero Mach number and Mach number equal to -0.5. 

Figure 4 shows comparisons at ka=20 of the reduction in sound pressure level versus polar angle 
averaged over azimuth with external sensors and internal sensors for and Mz=0 respectively using 18 
regularly spaced sources. 

Figure 5 shows the change in modal amplitude verses propagation angle following the minimisation of 
the sum of squared signal at ten equally spaced line arrays comprising ten elements, each forming beams 
in the directions between 60° and 90° in 5° increments, in the example in figure 4. 

Figure 6 shows comparisons at ka=20 of the reduction in sound pressure level versus polar angle 
averaged over azimuth with external sensors and internal sensors for M = -0.5 respectively using 18 
regularly spaced sources. 

Figure 7 shows the change in modal amplitude verses propagation angle following the minimisation of 
the sum of squared signal at ten equally spaced line arrays comprising ten elements, each forming beams 
in the directions between 30° and 60° in 5° increments, for the example of figure 5. 

Figure 8 shows the relationship between phase velocity, croup velocity and intake axial flow velocity in 
aduct. 

Figure 9 shows an unflanged hard walled duct containing a subsonic intake flow. 
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Noise from a.rcraft cause annoyance in populated areas. For an aircraft on approach and during takeoff 
these noise radiation angles are well away from the duct axis, towards the sidelines. Reference to the 
term "control bandwidth" means the band of angles (measured from 90° to the duct axis, over which the 
radiated sound is to be minimised. The relationship between transmitted and radiated sound fields has 
been determined, which allow a method of active control to be formulated in order to reduce the radiated 
sound power in a band of angles towards the sidelines. Matching the dimensions of the control region to 
the typical beamwidths of the principal far field radiation lobes produces reductions in the transmitted 
sound field in a continuous band of propagation angles. Modes whose main radiation lobe is in this band 
are nearest cut - off and are characterised by steep propagation angles relative to the duct axis. 

the invention uses e.g. one or more circumferential arrays of appropriately phased sensors at the duct 
wall that can observe the acoustic pressure associated with those propagation angles that are responsible 
for the field radiated towards the sidelines. 

Before describing the implementation of the invention, the fofrnulation of the relationship between 
reductions in the radiated far field and reductions in the internal, transmitted field are shown and terms 
are hereinafter defined. Most of the other terminology is familiar to the person skilled in the art. 

Figure la represents a circular, hard walled flanged duct ( I ) containing uniform, axial flow of Mach No. 
of M, (2). A ray-mode of acoustic pressure (sound) (2) transmitted along the duct and then radiated from 
the duct intake is shown and can be detected by a wall mounted line array, of appropriately phased 
(located) sensors (3). The figure shows the relationship between resultant wavenumber k m(1 normal to 
the local wavefront, the propagation angle 6 mn and axial wavenumber k OTn . The symbol k designates 
acoustic wavenumber k=2nf/c where f is the sound frequency and c is the speed of sound. The in-duct 
error" sensing principle proposed here is based on the mode angle Q mn which specifies the angle between 
the modal wavefront and the duct axis. More importantly the mode angle 9 m „ in the duct, for both 
flanged and unflanged ducts, is also coincident with the angle of the principal lobe of far field radiation 
providing there is zero flow external to the duct. Even when the flow speeds inside and outside the duct 
are different, a unique and monotonic relationship exists between the transmission and radiation angle. 
Much of the original mode-ray angle information pertaining to the transmitted sound field inside the duct 
is therefore preserved in the radiated sound field. The duct also contains one or more secondany sources 
(4) preferably, as angular arrays, the control of which is dependent upon the received signals of the 



sensors. 



The sensors are arranged in the figure as a series of annular rings. However the method is not limited to 
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such an arrangement, and may include any suitable spacing e.g. the sensors may not form rings but ma> 
be clustered closer together over a small sector of the duct wall. 

Initial studies revealed the performance of a single circumferential array of up to 20 wall mounted 
secondary sources whose strengths were determined in order to minimise the sound power radiated from 
the intake that passes across a hypothetical far field surface which subtends a band of radiation angles to 
the duct axis. Performance predictions were obtained over a frequency range of 0 < ka < 25 for a 6m 
long duct of 1.5m radius with intake flow corresponding to a Mach number KL = -0.5. The secondary 
source an-ay was located at a single axial location 2m from the face of the duct intake. Reduction in the 
amplitude of modes with the greatest propagation angle are responsible for sound power reductions in 
this desired band of radiation angles that are directed towards the sidelines. The different responses of 
the modal amplitudes found with and without flow can be attributed to an increased number of modes 
that can propagate in a mean flow and to a diminished range of propagation angles that follow compared 
to that no mean flow. 

It was deduced that increased number of modes can propagate in axial, mean flow compared with no 
flow, and this results in diminished range of propagation angles. By contrast, with flow, the amplitudes 
of the majority of modes that propagate close to the duct axis, such as the plane wave, are increased, 
resulting in an increase in sound pressure level at radiation angles in the forward directions. These are 
not a significant contributor to community annoyance due to the very long propagation distances to the 
ground. Initial modelling detennined that secondary sources are not required to reproduce sensitive 
phase changes which affect a number of modes. 

The following mathematically quantifies the relationship between the modal axial wavenumber and the 
angle between the modal wavefront and the duct axis in circular, hard walled ducts containing a uniform, 
axial flow. This relationship will be found in later sections to be central to the design of the in-duct 
sensor array. The mathematics also sheds light on the preferred design features of the invention. 

The acoustic pressure in a circular duct satisfies the convected form of the wave equation written below 
in cylindrical co-ordinates. 

(d : i a- aM - \ M c ., ay , 

1 ^ " 7^ ^ ^ ^t]p(^*. -0= : - A/ r —J ,(a.o.r.f) m 
where the coordinate system is defined in figure I. SI. is the Mach number of the flow and c is the 
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ambieni sound speed. This wave equation is defined such that A/ : < 0 at the due, intake. At sufficiently 
high ka. reflected sound at the duct termination is negligible. The complete solution to this equation for 
harmonicallv time varying sources in a circular flow duct, neglecting reflections, has the separable form 



(2) 



where and ^denote the modal amplitude and the axial wavenumber associated with modes 
propagatilg towards the duct exit and z, is the axial location of the source. In a hard walled duct the 
radial eigenvalues k mn equals; m /a , where; m „ denotes the n' h zero of, m „ and a is the duct radius. 

Substituting this solution back into the wave equation yields the following dispersion relat.onship 

where k is the free space wavenumber »/c. k^, actually represents a combined radial-circumferential 
wavenumber. The resultant wavenumber in the duct k mn is therefore equal to k - MJc^ . The angle 9^, . 
which specifies the angle between the modal wavefront and the duct axis, is calculated from cos8» = 
Arzm n /*m»thus 

„ ~ _ *z~ (4) 

Equation (4) can be re-arranged to express the axial wavenumber of the ( m ,f mode in terms of the 
propagation angle as 

, _ -fecose,,,, ( 5a) 
z """l + M.cose„,„ 



k ~ k (5b) 

- l+A/ : cose MM 

A geometric interpretation of these wavenumbers and their relationship to the modal propagation angle is 
illustrated in figure lb. A surprising aspect of equation (5b) is its independence of the duct radius a. Thts 
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is not the case for the ray mode angles in the radial and circumferential directions which vary quite 
strongly with a. 



The external, far field acoustic pressure due to the im.n/ h mode from a flanged duct may be written in 
the form 



. exp {/(aw + m§ - kR)\ 
Puui = D m „{k$) HVV R }J (6) 



where D mft denotes the directivity function of the fiw.ii/* mode and R is the distance from the duct face to 
the observer. 



n jk(j sin*)- J UtosmQ)) 



Expressing 9^ in terms of the axial wavenumber with the aid of the dispersion relation of equation (3) 
for M : = 0 gives 

6 p^cos- {k u Jk) (8) 



This describes the essential monotonic and unique relationship between the axial propagation angle Q m „ 
of a mode in the duct. 

A wall mounted phased line array for the detection of modes by modal angle 8 mn . Modal amplitude 
reductions that result from reducing the sound power radiated in a band of angles towards the sidelines 
has been shown to bear a definite and causal relationship to reductions in the in-duct sound field 
transmitted obliquely to the duct axis. The objective is therefore to design a wall mounted sensor array 
comprising of a relatively small number of discrete sensors that has sufficiently good directivity to detect 
this change in the transmitted sound field. Since modes can only be controlled if they can be observed, 
the ideal receiver response is plotted in figure 4 and is a step function which detects only the signals 
arriving at large incidence angles to the array while rejecting signals transmitted at angles close to the 
duct axis. Also plotted in figure 4 is the typical variation of axial wavenumber with propagation angle 
from equation (5a). Figure 4 therefore demonstrates that the ideal receiver characteristics is a high pass 
filter ot propagation angle which, by virtue of equation (5). is also a low pass filter of axial wavenumber. 
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It shows a rypical variation of axial wavenumber with propagation angle and the ideal receiver response. 

A simple sensor array whose directivity characteristics approximates to the ideal step function response 
illustrated in figure 4 will now be described. The relationship between axial wavenumber and 
propagation angle given by equation (5a) can be used to express the acoustic pressure at some 
circumferential location «|» at the duct wall r = o such that 



The pressure signals at an axial line array of sensors at the duct wall due to the internal transmitted sound 
field is therefore indistinguishable from a series of plane waves arriving at the modal propagation angles. 
To be able to discriminate all possible modal arrival angles (or more precisely spatial frequency) without 
ambiguity, the minimum sensor separation distance Dz is required to be. as a consequence of the Nyquist 
sampling theorem, one half the wavelength of the highest axial spatial frequency in the flow. According 
to equation (5a), this is the plane wave mode (corresponding to 8 M = 0) transmitted at the highest 
frequency of interest in the highest intake flow speed of interest M^. The highest frequency in 
general corresponds to the highest harmonic s maj of the blade passing frequency of interest = s^Clb. 
where O is the shaft rotational frequency, and b is the number of rotor blades. From equation (5a), 

*-**«.0 + "i-) (10) 

where - Kmm = c/f^ which is the shortest wavelength in the radiation field. The frequency, /», is known 
as the design frequency of the array. If z, denotes the axial position of the first sensor in a wall mounted 
line array comprising I elements separated by a distance te. the f sensor is required to have the axial 
location z t given by, 

r,=z (+ (/-l)Az HI) 

the acoustic pressure ptaAzj at the t" sensor can therefore be written in the form 
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where 6^ is simply a phase term that is constant across the array, thus 



1 + iW.cose^ 



and y mn is the relative phase angle between adjacent sensors 



*{n/„X* (l4) 

l+w.cose^ 



In order to prevent ambiguity in the measured arrival angle caused by aliasing, therefore, a 
fundamentally important condition is that / < /„. .Equation (12) effectively specifies the complex 
weights wifaaJ of a simple 'delay and add* line array beam former. In order to preferentially amplify the 
acoustic pressure signal arriving at angle Q 0 to the line array,, the anay elements are simply required to 
delay the signals at each sensor by an appropriate amount Vofia) which upon addition, causes the signals 
at each sensor to be summed perfectly in-phase. The beam steer angle are made such that they are made 
to scan the angles 6 from 90° to 90° to 90° -A8 in some appropriate incremental angle. By inspection of 
equation ( 12), therefore, the f h element in the line array is required to have the phase specified by 

w,(e 0 )=expfy/u/ 0 (e 0 )} <is> 

where 9 0 is the beam steer angle 



^(///^X^A/ 2w Jco $ e 0 

• 1 + jW.cosG 0 



This formula is a fundamental result which enables the invention to be implemented and allows a method 
for devising a directional wall mounted receiver simply by locating the axial line array of sensors with a 
maximum separation distance equal to Ar = X mm (l +. W^) and by introducing the relative time 
delays between the sensors specified by equations (19) and (20). It is generally applicable to all circular 
ducts, irrespective of radius and depends only on the Mach number of the free stream flow and the 
frequency of the fan tone to be controlled. Both these parameters can be readily determined. 
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Sensor line array directivity (beam Steer angles) 

The following describes relationships between the sensor (receiver) line array directivity characteristics 
and its relationship to mode detection. 

The directivity characteristics of the sensor (receiver) line array can be described by the normalised 
directivity function rf/B^- This function specifies the array response at some angle 9 when the main 
beam is steered at an angle 0^, and can be determined from 



4e/e 0 )=-5-texpfy/(w 0 )} (") 



and is defined such that dtoSnt = I. Equation (21) is a geometric series that can be summed over L 
terms to give 



#n,U *y -"'11} 

v °' I sin((vv =v.)/2) 



(18) 



In terms of the propagation angle, v - u/„, mav be written as 

v I c os6 cos9 0 N | 

V " V , = * (/// m „Xl + \ j + M cos9 " , + M : cos9 0 J ( } 

By way of example, comparisons of the directivity functions evaluated at the design frequency J m<a at 
Mach numbers of 0 and -0.5 for an amy comprising ten elements steered at 45° to the array axis is 
presented in figure 5. Figure 3 shows the directivity function of a ten element line array at the design 
frequency, steered at 45° for zero Mach number (solid line) and with a Mach number equal to -0.5 
(dashed line). For this particular receiver array the lOdB beamwidth is about 20°. The presence of flow 
with speed equal to M. = -0.5. which is typical for an aircraft on approach, appears to cause no 
appreciable change in directivity characteristics apart from a slight narrow.ng of the main beam and a 
reduction in the number of side lobes. The important difference is that to achieve roughly ihe same 
beamwidth in this now speed, a ten element array at the design frequency has a shorter length of just 
2.25/.„„„ which is half the array length necessary in a duct without flow according to equation (14). For 
this reason it would appear, therefore, that the presence of intake now. having the property of contracting 
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ihe spaiial frequencies by virtue of an effectively reduced sound speed, is advantageous. For example, at 
the comparatively low frequency of 500Hz. [tea = 13 for a = 1.5m) the optimal length of a ten element 
array at the duct intake at the design frequency is about 1.5m. A reasonably long array with good 
directivity can therefore be fitted, quite easily, within modern high by - pass ratio engines. At more 
realistic blade passing frequencies, that are typically greater than 500Hz. the array lengh is even shorter. 

The output b(Qn.a,$) of a receiver array located at an azimuthal angle 6 at the duct wall r = a and steered 
at an angle Q n is obtained from the following operation on the discrete wall pressure measurements 

^..^)-Z w ^oWM^/) (2°) 
/-i 

where the acoustic pressure at each sensor is the sum of ray - modes according to equation ( 12) so that 

^..fl.O-Zp-M.^)«v{> 5 -}Z ex p(/ / ^--*o)} (20 

The summation of terms over / can be evaluated exactly to give 

b(Q o , a, <j. ) = £]T — ^ r--f p m „ (a.<t>, z, )exp j5 M (22) 

which is precisely the directivity function Ld(Q m ^Qn) of the receiver line array deduced in equation (22). 
The receiver output can therefore be written as 

b(Q 0 ,a.^)=LY J 4^^o)pU^- = ,>MjOn U ,} (23) 

HJ.M 

The effect of implementing this receiver array is to weight the modal contributions to the receiver by a 
factor equal to the array's directivity function df$„ t JQ n t evaluated at the modal arrival angle 8„ m . 
Steering of the array s main beam in the direction of the mode angles closest to cut-off will therefore 
amplify the acoustic pressure propagating with those angles highlighted in figures 2 and 3 as being 
directly responsible for the reductions in the important band of radiation angles, i.e. those towards the 
sidelines. The transmined sound field whose propagation angles are diffracted outside the control 
reeion. i.e.. close to the duct axis, will be partially rejected at the receiver by an amount depending on the 
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degree of side lobe suppression compared to the main beam according to equation (23). The ability of the 
receiver array to discriminate between different modal propagation angles will depend on the angular 
bandwidth of the transmitted sound field in which pressure reductions arise, in relation to the beamwidth 
of the array. Greater rejection of the signals due to the unwanted arrival angles will be achieved when 
the ray modes are well separated in propagation angles. 



Implementation the inventi o n for far field active control USinr in-dltCt fine arT3V gfTPf SenSOfS - 

The above describes relationships between the in-duct pressure field and how they effect the acoustic 
pressure in the far field. The above analysis has shown how best to glean in-duct sensor data to estimate 
far field effects. This enables the skilled man to optimally design the sensor array including beam steer 
angles and the use of general parameter such as duct flow speed to estimate far field radiated pressure. 
The important point here is the ability to quantify the effect of changes in the in-duct field to changes in 
the far field. From initial calculations sound power reductions at radiation angles towards the sidelines is 
accompanied by well defined changes to the transmitted sound field and that the change in the angular 
variation of the transmitted field was detectable by a line array receiver at the duct wall comprising 
relatively few sensors. 

In this section the use of these sources to implement the invention is described whereby control of 
secondary sources (loudspeakers) arranged in the duct are suitably operated from data of the sensors to 
optimally minimise sound at a particular fan tone frequency. The secondary sources will be driven to 
minimise the sum of squared signals each signal being produced by steering a beam formed by a number 
of independent axial sensors line array located around the duct wall: this brings about a similar 
modification to the transmitted field as that produced by conventional external far field error sensors, in 
order to procure similar reductions in the radiated field. This is done by using the algorithms set out 
below. 

If K line arrays each produces / signals by steering of the beams at / angles, a suitable cost function J is 
given by 

.-I *.i 

where W9 0 ,.u.<h> denotes the complex signal produced after steering a beam at an angle 6„, by a receiver 
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array located at the circumferential angle ^ around the duct wall and is computed from 

/. 

*( 9 o,^-*i)=£^^o,>Mi--i) (25) 

From hereon the dimensions of the wall mounted sensor array will be denoted by (K.L) so that the total 
number of error sensors is K x L. The function J in equation (24) call be expanded to produce a 
quadratic function of the secondary source strengths. The vector of optimal secondany source strengths 
that uniquely and globally minimises J can be deduced by standard mathematical methods. Note that it 
is desirable to have as many sensors as possible in the axial direction to provide good receiver 
directivity, and as many line arrays around the duct wall as possible to ensure that J is minimised. 

Figure 4 shows tests of the in-duct receiver array's ability to control the radiated sound towards the 
sidelines is from a duct without flow. A 10 x 10 sensor amy is used comprising ten line array receivers 
equally spaced around the duct wall, each consisting of ten elements. The beams at each of the receivers 
are steered in the range of angles from 55° to 90° from the duct axis in increments of 5°. Eighteen 
secondary sources are driven to minimise the sum of squared signals produced by the ten independent 
receivers according to equation (24). A comparison between the radiated far field sound pressure level 
reductions, obtained by computer simulation versus polar angle produced by using the internal and 
external sensors is shown. These results represent the average reduction over twenty azimuthal angles. 
The solid curve is the result of minimising the sound power radiated into a band of angles from 55° to 
90° from the duct axis using a dense grid of external error sensors in the control region that afford 
perfect observability of the radiated field. As a result of the measurements from the sensors located 
inside the duct, loudspeaker or sources can be driven to minimise the noise in the far field. Figure 5 
shows the change in modal amplitude verses propagation angle following the minimisation of the sum of 
squared signal at ten equally spaced line arrays comprising ten elements, each forming beams in the 
directions between 60° and 90° in 5° increments. The agreement between the two curves is extremely 
good. Sound power reduction in the control region using the in-duct receiver array is !7.9dB which 
compares very well with the theoretical maximum of 2l.6dB obtained by the external sensors. Both in- 
duct and external error sensing strategies produce the necessary, and very similar transformations in the 
transmitted sound field. 

Increasing the dimensions of the array to (20.20). requiring a total of 400 microphones, affords a further 
increase in sound power reduction to 20.6dB. Although the number of sensors is unrealistically large in 
this case, this simulation serves to validate further the principle underlying the proposed control 
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technique. Howev e r. increasing the number of sensors indefinitely does no. produce reductions 
arbitrarily close to the theoretical maximum since one is ultimately limited by the optimal phase 
relationships berween the modes. The proposed sensing technique does not allow this manipulation of 
the modes. The significance of the phase in the control mechan.sm. while comparatively unimportant at 
high ka. is important when the number of propagat.ng modes is small, i.e.. at low ka. Nevertheless, 
substantial pressure reductions are still achievable at these low frequencies. 

The same performance comparison of the change in sound pressure level versus polar angle berween 
using the internal and external error sensors was repeated at ka = 1 5 with a flow speed of hU = -0.5 and is 
presented in figure 6. A ring of fifteen secondary sources were driven to minimise J comprising signals 
produced by the (.0.10) sensor array forming beams steered at 30° to 60° to the duct axis in 5° 
increments. These beam steer angles differ from the previous no flow example to take account of the 
modified propaaation angles due to the flow. As before, the two results agree to an extremely good 
degree. The sound power reduction in the control bandwidth is 10.5dB which is just 3dB below 
theoretical maximum reduction of 13.4dB obtained using far field sensors. However, the effect on the 
• transmitted field produced by the two approaches is slightly different, although the broad mechanism of 
control which consists of reducing the modal amplitudes closest to cut-off remains the same. F.gure 7 
shows the change in modal amplitude verses propagation angle following the minimisation of the sum of 
squared signal at ten equally spaced line arrays comprising ten elements, each forming beams in the 
directions benveen 30° and 60° in 5° increments, for the case in figure 5. 

A large number of sensors is required to construct the array. About one hundred is anticipated to be 
necessarv, although many more would of course be des.rab.e. However the large number of sensors 
required does not translate to a correspondingly high processing bandwidth. The reason for this ,s the 
independence of the line arrays which is fundamental to ensuring robustness of the technique. One ,s not 
required to measure transfer functions between sensors on different line arrays. Each receiver line array 
could therefore be allocated its own dedicated processor for forming the beams whose output could then 
be input to a main processor for real time adaptation of the secondary sources. 

Another important constraint on the receiver array arises from the importance of sensing only the in-duct 
propagate field, which contains all the information contained in the far field radiation, and not the 
evanescent field close to the sources which transports very little energy to the radiated far field. The 
sensor Cement nearest the sources must therefore be separated by several acoustic wavelengths in order 
to avoid contamination of the measured signals by the non-propagating field. 
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A further issue relates to which beam steer angles gives best reductions in far field noise level. The 
relationship between the propagation and radiation angles in real turbo fan engines svill certainly be 
much more complicated than that suggested here: ray mode angles are complicated by complex 
geometry of the nacelle, shear velocity' and temperature profiles present across the duct. However in all 
cases from elementary acoustics a unique, monotonic correspondence exits between the angle of the ray 
mode at the duct and the radiated angle of the far field peak pressure maximum, allowing the in-duct 
sensing technique of the invention to be successful and widely applicable. 

By minimising the sum of squared signals produced by a number of identical, independent line arrays 
equally spaced around the duct wall, each forming beams at the appropriate angles, a similar 
modification to the transmitted sound field (radiated far field) is obtained by the invention. Using the 
internal sensor array gave reductions in the radiated sound power tosvards the sidelines which was within 
a few decibels of the theoretical maximum reduction obtained given perfect observability of the radiated 
field. Significantly, these sound power reductions were achieved without the knowledge of transfer 
functions between sensors on different line arrays; the technique is therefore likely to be stable and 
robust by virtue of its simplicity. 

The beam steer angles that afford the best reduction in radiated pressure towards the sidelines are 
therefore very difficult to predict in real turbofans. However, it is envisaged that in practice the best 
combination of beam steer angles, number of secondary sources and number of line arTays etc. which 
afford the greatest reduction in noise will be determined from the results of a number of systematic fly- 
by tests. 

The sensor array should be preferably located flush to the duct walls in order not to interfere with the 
passage of flow through the engine. 

Relationship between in-duct propagation angle and the far field radiation angle 

The use of the sensor array is fundamentally dependent on the existence of a unique and simple 
relationship between the transmitted sound field and the radiated far field. The cost functions to be 
minimized depend on whether total sound power or sound pressure towards the sidelenes is to be 
reduced. In order to enable these cost functions, which are given hereinafter and hereinbefore, to be 
implemented, basic relationships and definitions are given. The sensor array described here is designed 
to detect the modes based on their different axial propagation angles 0 mp . By simple geometry this angle 
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is also given by: 



c 



0 = — 



(26) 



where c m denotes the modal phase speed and may be regarded as a vector normal to the modal 
wavefront. A more fundamental variable is the angle with which acoustic energy is transmitted along the 
duct and this is related to the axial group velocity c gmu . where 



. = JUL. (27) 

Perfoming the differentiation of the dispersion relation of equation (3) yields the following relationship 
between the axial phase and group velocities 



c« = c.„ +c 0 M : 



(28) 



The angle with which acoustic energy is transmitted along the duct is identical to the angle of the mode 
peak pressure far field radiation lobe 8, m , when the flow speed inside and outside the duct are equal. 
Thus, 



o c ,..» +c o M : (29) 



where c R is the resultant sound speed in the direction of sound power propagation. This result is readily 
derived since the sound speed in the radial and circumferential direction are unchanged by the presence 
of fiow. A sketch indicating the relationship between the phase and group velocities is provided by 
figure 8 for the case when the flow speed is the same everywhere. 

The angle Q? mu of the modal peak pressure radiation 

- _ k ^ ~ M ': (30) 
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cost) 



(31) 



which allows comparison with the expression for 9 mM in terms of the cut-off ratio 




(32) 



Thus, the in-duct axial propagation angle 8 mM which can be detected using the in-duct sensor array, and 
the radiated peak far field pressure angle 8/> mM , are connected by the above two equations with cut-off 
ratio as the independent parameter. 

The angle of the modal peak far field pressure radiation lobe is closely related to the axial propagation 
angle. The deviation between the two angles increases with increasing Mach number. The range of axial 
propagation angles becomes smaller with increasing intake Mach number whereas the range of the 
principal radiation lobe angles remains distributed between 0° for the plane wave mode and 90° at cut- 
off, irrespective of flow speed. 

Reducing transmitted and radiated sound power 

The active minimisation of the total transmitted sound power was found to be particularly effective in 
reducing levels at the fundamental of the blade passing frequency. However, sound power was not found 
to be an appropriate cost function when reducing sideline radiation is the main objective, although there 
may be occasions when reducing the sound power may be desirable. For uncoupled modes the sound 
power transmitted in the flow duct could be written as: 



ft' = y v \b |\ y 



>0 



(33) 



where •/„,„ is the real pan of an effective modal admittance in the flow calculated from 
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where (3 WM is the modal admittance 

The real part of the effective modal admittance in the flow appears to be closely correlated to the modal 
propagation angle 9 mM , and tends to zero as the mode approaches cut-off. ie.. y mM -> 0 as 8 mM -* 9. wat . 
No sound power is therefore transmitted at cut-off. A good approximation to the effective modal 
admittance y (9) plotted in figure 3 versus propagation angle which incorporates this cut-off condition is 
described by the simple function 

Y(0)=cos(i7t0/e 2m „) 05) 

where Q 2tmu = cos\-M z ). A cost function approximately proportional to the transmitted sound power is 
therefore given by 



A-itf$.w>i.«.*.f (36) 

A-l i»l 

where 9, spans the full range propagation angles such that 8, = 0° and 9/ = Q snax in some appropriate 
incremental angle. For completeness, recall that 



*#..«.♦)- Z*i$o>M.*i) (37) 

/-i 

and 

w,(e 0 )=exp{;7 M /(e 0 )} (38) 



V ^ l+A/ 2 cos9 



(39) 



An imponant feature of the cost function of equation (36) is the summation over axial arrays at different 
azimuthal positions around the duct wall. A cost function based on a single axial line array would, in the 
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general case, be minimised by rotating various spinning modes to produce destructive interference at the 
azimuthal location coinciding with the location of the sensor array. Minimising the sum of square 
outputs from several axial arrays prevents this from occurring and ensures the correct control mechanism 
by reducing the appropriate modal amplitudes. The number arrays A' should be made equal to the 
number of significant circumferential modal order present in the duct. However, when only a single 
modal order m is present, is often the case at the comparatively low frequency corresponding to I BPF. a 
single axial sensor array is entirely sufficient and A' = I . 

Reducing sound pressure levels towards the sidelines 

In order to reduce by active means the sideline radiation it is necessary to target these radiation angles 
specifically rather than minimise a global quantity such as sound power. Furthermore, azimuthal 
directiviry of fan noise is generally weak by virtue of the very small number of circumferential modal 
orders which are able to cut on. Implementing active control in a narrow band of polar radiation angles, 
but extending the control region to include all azimuthal angles, incurs no appreciable performance 
penalty. Hereinafter is described a method of how the in-duct error sensor array provides sound pressure 
level reductions over an ax i-sym metric control surface in a band of polar radiation angles. A0. 

The control mechanism underlying the reduction of sideline radiation is made clear by the relationship 
between the in-duct angles and those in the far field. Tie modes closest to cut-off must be attenuated in 
order to reduce the sideline radiation. A sketch of the duct, the control surface and axial sensor array 
configured to control sideline radiation is presented in figure 5 showing an unflanged circular hard 
walled duct containing a mean subsonic intake flow. Enclosing the duct exit is a sector of a sphere of 
width A9 across which the sound power is to be minimised by a ring of secondary sources. A single 
axial sensor line array at the duct wall detects the transmitted sound field from which the field radiated 
towards the sidelines can be inferred. 

A suitable weighting function on \b(B h a$ k )\ 2 which accounts for the importance to sideline radiation of 
near cut-off modes is the exponential function exp(-n(6 :mai -8)}. where rj is an arbitrar> constant which 
specifies the relative weighting assigned to the different modes according to propagation angle. A cost 
function for reducing sideline radiation incorporating this weighting function is expressed 

J *=t "£«PH^«.-e,)|A#,.a.4 t J : - (40) 
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The range of 'look* angles 9 is now taken over the range of propagation angles close to 9. mai that are most 
responsible for the sideline radiation. Note, that the summation over different azimuthal axial sensor 
arrays K is included for generality in order to allow for the presence of a number of different 
circumferential modal orders m. 
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Claims 



I A duct for fluid flow having means for active control of sound radiated therefrom, said duct 
comprising sound sensors located on the inner surface of said duct and grouped together in one or more 
planes transverse with respect to the duct axis, and at least one secondary source whose operation is a 
function of sound received at said sound sensors characterised in that the axial spacing of said transverse 
planes is not more than 0.5UO + w *ere is the wavelength corresponding to the radiated 

tone frequency of interest and M zmax is the maximum Mach number of the free stream flow in the duct. 

2. A duct as claimed in any of the above claims characterised in thai said secondary sources are 
located on the inner surface of said duct and grouped together in one or more planes transverse with 
respect to the duct. 

3. A duct as claimed in any of the above claims characterised in that the secondary sources are 
controlled so as to minimise the received sensor signals in a pre-set band of angles to the duct. 

4. A duct as claimed in any of the above claims characterised in that the sound sensors have a step 
response whereby only sound propagated within a pre-set angle to the duct axis is effective sensed and 
used to control the operation of said secondary sources. 

5. A duct as claimed in any of the above claims characterised in that the secondary sources are 
operatively controlled so as to minimise the sum of squared signals received at the sensors. 

6. A duct as claimed in any of the above claims wherein the secondary sources are controlled as a 
function of the sensors so as to minimise the cost function 



where b(Q fl ,M$k) denotes the complex signal produced after steering a beam at an angle 8*, by a receiver 
array located at the circumferential angle $ k around the duct wall and is computed from 



(41) 



^e 0 ,MO=Z^(eo>(Mv--,) 



(42) 
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uv(9 0 )=exp{-/7\j/ 0 1 (43.) 



where y 0 (9") is the beam steer angle 



Vo=— — : — 77 (44) 

1 + m z cos0 o 



7 A duct as claimed in any of claims I to 5 wherein the secondary sources are controlled as a 
function ofthe sound sources so as to minimise the cost function 

A>l i*l 

where 

y(e)=cos(lne/9 2miJl ) (46) 

e. m «=cos-'(-M z ) (47) 

8. A duct as claimed in any of claims I to 5 wherein the secondary sources are controlled as a 
function of the sound sources so as to minimise the cost function 



/'aft- 1 Zexp{-i!(e 2in „-e,)56(e„a,^]| 2 (48) 
i«io»e f >e / ^. 



r\ is an arbitrary constant 

9. A method for the active control of sound radiated from a fluid flow duct comprising: 

a) sensing sound from an array of sensors located on the inside surface of the duct and 

b) controlling an array of secondary sources located on the inside surface of the duct so as to minimise 
sound radiated in the far field in a pre-set band of angles to the duct axis, characterised in that only sound 
propagated within a pre-set angular interval to the duct axis is effectively sensed and used to control the 
operation of said secondary sources. 
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10. A method as claimed in claim 9 characterised in that step (b) comprises controlling an array of 
secondary* sources so as to minimise the sum of squared signals received at the sensors within a pre-set 
angular interval. 

11. A method as claimed in claims 9 or 10 characterised in that the secondary sources are controlled 
so as to minimise the cost function: 

where i>(8o,.tf.<t>*) denotes the complex signal produced after steering a beam at an angle 6 0l by a receiver 
array located at the circumferential angle $ k around the duct wall and is computed from 

w,(e o )=exp{-,7y 0 } (SI) 
where y^Qo) is the beam steer angle 



0 1+M 2 cos9 0 



12. A method as claimed in claims 9 or 10 wherein the secondary sources are controlled to minimise 
the cost function: 

*-l i.l 

where 

•;(e)=cos(tne/e 2in j m 
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8,.. -«*-(-*,) (55, 

13. A method as claimed in claims 9 or 10 wherein the secondary sources are controlled to 

minimise the cost function: 

r| is an arbitrary constant 
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Fig. 1b. 




(b) 
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Fig.2. 
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Fig.4. 
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Fig.6. 
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